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Coherent state functional integrals for the minisuperspace models of quantum cosmology are 
studied. By the well-established canonical theories, the transition amplitudes in the path-integral 
representations of Wheeler- De Witt quantum cosmology and loop quantum cosmology can be for- 
mulated through group averaging. The effective action and Hamiltonian with higher-order quantum 
corrections are thus obtained in both models within the scheme of Gaussian coherent states. It turns 
out that for a non-symmetric Ifamiltonian constraint operator, the Moyal (star) -product emerges 
naturally in the effective Hamiltonian. This reveals the intrinsic relation among coherent state 
functional integral, effective theory and star-product. Moreover, both the resulted effective theories 
(N- imply a possible quantum cosmological effect in large scale limit under certain condition. 
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I. INTRODUCTION 



One of the most fundamental issues in modern physics is quantum gravity. Among various approaches to quantum 
gravity, the viewpoint of nonperturbative quantization has been received increased attention. The first significant 
attempt of this kind was presented in 1960s [1]. Based on the ADM Hamiltonian formalism of general relativity (GR) 
and the Dirac's generalized Hamiltonian quantization theory, which is in principle applicable to constrained system. 
Wheeler and DeWitt proposed a wave functional description of gravitational field. In this quantum geometrodynamical 
approach, the whole universe is described by a wave function which is defined on superspace - the space of all 3-metrics 
and matter field configurations. The dynamics is encoded in the second-order differential Wheeler-DeWitt (WDW) 
' equation. Despite the elegant form of WDW theory, it encountered a number of fundamental problems, such as the 
^ . physical meaning of quantum constraint equation, ordering of operators and problem of time [2] . To overcome these 
obstacles, especially to impose boundary conditions on WDW equation to find solutions, there were some attempts 
appealing to the path integral formalism. Hartle and Hawking suggested an Euclidean path integral representation 
ly-j ■ [3], which gave the "no-boundary" proposal for the wave function. However, this approach can not be evaluated 
• ] exactly in general cases as the integral was usually badly divergent. To alleviate the difficulty, Halliwell derived a 

■ Lorenzian path integral formalism for minisuperspace models [4], which revealed the relationship between the choice 
I ■ of measure in the path integral and the operator ordering in the WDW equation. 

, In the last two decades, an alternative background independent approach developing rapidly is loop quantum 

■ gravity (LQG) [5-8]. The starting point of LQG is the Hamiltonian connection dynamics of GR rather than the ADM 
. ^ , formalism. In this framework, GR looks like a gauge field theory with SU{2) as its internal gauge group. By taking the 

■ holonomy of su(2)-connection A\ and flux of densitized triad i?^ as basic variables, the quantum kinematical framework 
5h [ of LQG has been rigorously constructed, and the Hamiltonian constraint operator can also be well defined to represent 

quantum dynamics. Moreover, a few physically significant results, especially the resolution of big bang singularity, 
have been obtained in the minisuperspace models of loop quantum cosmology (LQG) [9, 10]. The quantum bounce 
replacement of big bang and its properties are being studied from different prescriptions of LQG [11-16]. Effective 
equations were also derived in isotropic models [19-21], which predict evolution of universe with quantum corrections 
and shed new lights on the singularity resolution. We refer to [17, 18] for overviews on the recent progress in LQG. 
Besides the canonical formalism, the so-called spin foam models were proposed as the path integral formalism of LQG 
[5]. However, whether the two approaches are equivalent to each other is a longstanding open question. Thanks to the 
development of LQG, we have much simple models to address this question. As symmetry-reduced models, there are 
only finite number of degrees of freedom in LQG. Following the conventional method in quantum mechanics, one can 
find the path integral formalism of LQG starting with the canonical formalism. This approach is being implemented 
by a series of papers [22, 23] with the scheme of simplified LQG [24]. In these papers, the transition amplitude 
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were castcd into path integrals by using either de-parameterized Hamiltonian or the group averaging method. Here 
one employed the complete basis of eigen-states of the volume operator to formulate a path integral with somehow 
descrete-steps, which inherited certain properties of spin foams [25]. Moreover, the first-order effective action for 
the path integral was derived by this approach [22, 23], which implied the origin of singularity resolution of LQC in 
the path integral representation. In canonical LQC, the effective Hamiltonian constraint with higher-order quantum 
corrections could even be obtained by the semiclassical analysis using coherent states of Gaussian type, which implied 
a possible effect of quantum gravity on large scale cosmology [11, 21]. Moreover, a systematic approach to the effective 
theories of quantum cosmology in the canonical framework is being developed [26, 27]. It is thus interesting to see 
whether the effective Hamiltonian in the canonical theories can be confirmed by some path integral representation. 
Since the higher-order corrections of the Hamiltonian come from the quantum fluctuations, a natural attempt to 
achieve them is to employ coherent state path integral [28] . 

In both WDW quantum cosmology and LQC, by coupling with a massless scalar field, the Hamiltonian constraint 
equations can be reformulated as Klein-Gordon-like equations, where the corresponding gravitational Hamiltonian 
operators, as multiplications of several self-adjoint operators, are non-symmetric in the kinematical Hilbert space 
[10, 24]. While this treatment is essential in order to obtain the physical states satisfying the constraint equation, it 
also provides elegant physical models to examine the so-called Moyal *-product in quantum mechanics. At the very 
beginning, Moyal proposed the ^-product in order to clarify the role of statistical concepts in quantum mechanical 
system [29]. Then this idea were generalized to many situations including quantum spacetime itself. It provided the 
core concept of non- commutative geometry, one of promising and interesting new tools in the study of quantum field 
theory and quantum gravity [30]. In canonical quantum theories, the * -product can also be understood by coherent 
state approach [31]. In fact, on a coherent state, the expectation value of the multiplication of two non-commutative 
self-adjoint operators equals to the * -product of the expectation values of the two operators. Thus it is also possible 
and desirable to derive the * -product in coherent state functional integral approach. The above models of quantum 
cosmology provide a good arena to launch the desired investigation. Hence we will study the coherent state functional 
integrals in spatially flat isotropic FRW cosmology coupled with a massless scalar field <j> in both WDW and loop 
quantization approaches. It should be noted that we will restrict our scheme to Gaussian coherent states for the path 
integral representation. 



II. CANONICAL FRAMEWORKS OF QUANTUM COSMOLOGY 

The Hilbert-Einstein action for this model is given by: 

S = J d^xV^R +\j (2.1) 

In this spatially flat model, we flx a space-like sub-manifold S, which is topologically and equipped with Cartesian 

coordinates = 1, 2, 3), and a fiducial flat metric °qab given by "qabdx'^dx^ = dx\ + dx\ + dxi^. The physical 3-metric 
qab is then determined by a scale factor a satisfying ciah = a^"qab- It is convenient to introduce an elementary cell V 
and restrict all integrations to this cell since the spatial slice is non-compact. The volume of V with respect to °qab 
is denoted as V„ and the physical volume is V — a^Vo- Then the geometrical pair {a,Pa) can be used as canonical 
variables, where the conjugate momentum satisfies Pa oc aa. 

In order to study the WDW cosmology and LQC on the same footing, we employ the new canonical variables 
(A^, Ef) in both theories. Due to the homogeneity and isotropy, we can fix a set of orthonormal cotriad and triad 
, °ef) compatible with "qab and adapted to V. Then the cotriad oj^ which are orthonormal with respect to physical 
metric qab can be written as cj* = x<i°'^a, where x = 1 if has the same orientation as the fiducial °a;^ and % = — 1 
if the orientation is opposite. The basic canonical variables take the simple form [32] 

= cVo^^'^'^^l = pV^F-(2/3)Oe«. (2.2) 

The dynamical variables are thus reduced to (c,p) with the Poisson bracket: {c,p} = 87767/3, where 7 is the Barbero- 
Immirzi parameter. Following the jl-schem of "improved dynamics" [10]. the regulator p, used in holonomies is given 
hy p, = \J A/]p], where A = A^/Sn^lp is a minimum nonzciro eigenvalue of the area operator [33]. In order to do the 
semiclassical analysis, it is convenient to introduce new dimensionless conjugate variables [11, 24]: 

4- 3^G(7r7)i£p Pa ^ _ sgn{p)\p\i ^ xa^Vp 3) 
2 Wo a2' ' 27r-f£lVA A ■ (n-f)! i^' 

with the Poisson bracket {v, 6} = i , where the Planck length £p is given by ~ Gh. From the matter part of 
action (2.1), we can get the momentum of </> as = ° ^""^ and the Poisson bracket: {(j),p<j,} = 1- The kinematical 
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Hilbcrt space of the quantum theory is supposed to be a tensor product of the gravitational and matter parts. In 
WDW quantum cosmology, one employed the standard Schrodinger representation for both matter and gravity to 
construct Hilbcrt space 7^^,?^- The (generalized) orthonormal basis in T^^irP^ given by |w) (S) |^) (or denoted as 
\v,(j))) with the inner product: {v' , (/)' \v , (j)) = S{v' ,v)5{<j)' The fundamental operators act on a quantum state 
4'(w, <p) e "Hkii?^ 3,s in standard quantum mechanics. To obtain the physical states, one has to solve the quantum 
Hamiltonian constraint equation [24]: 

-CwDW*(^^,0) = (^-^+ewDW^'i>{v,(p) = {dl-127rG-vd^{vd^,))^{v,cl)) = 0. (2.4) 

However, in LQC, while the Schrodinger representation was still used for the matter, gravity was quantized by the 
polymer-like representation [32]. Thus quantum states in the gravitational Hilbcrt space of LQC are functions ^{v) 
of V with support on a countable number of points and with finite norm ||^!.||^ := [34]. Hence the inner 

product is defined by a Kronecker delta {v'\v) = The basic operators act on a quantum state ^{v,(j)) in the 

kinematical Hilbert space 'hC^^ as: 

t)*(w, (j)) = v^{v, (j)), e*^^{v, 4i) = *.(« + 1, (j)). (2.5) 
The quantum Hamiltonian constraint equation becomes 

-^QC • = + ®LQC j <t>) = 0, (2.6) 

where Blqc is a positive second-order difference operator defined by a simplified scheme as [24] : 

6lqc • Mv, (t>) = [{v + 2)Mv + ^,4>)- 2vMv, 4>) + {v - 2)^{v - 4, 4,)] . (2.7) 

Solutions to the constraint equations and their physical inner products can be obtained through the group averaging 
procedure. It is demonstrated in [22] that in a timeless framework the most important entity is the amplitude 

AaC 



A{vf,(pf;v^,(pi) -.^ da{vf,(j)f\e"^^\vi,(l)i), (2.8) 
t/ — 00 

which contains all the dynamical information. 

We are going to concern about coherent state functional integrals. The {generalized) coherent state of the matter 
part is labeled by a complex variable Zg := -^^{(l>o + ■%<^'^P<t>„) and defined by 

l^-^J := /" e-^^'^l^e*f*<'('^-'^°)|^), (2.9) 

J — CO 

which is the eigenstate of the annihilation operator z = ■^^{4> + jiP4>'^'^), where a describes the width of the wave- 
packet or quantum fluctuation. It satisfies the key properties of a coherent state, namely, saturation of Heisenberg's 
uncertainty relation, resolution of identity and peakness property. Similarly, we can introduce rjo '■= ;^^(''^o + iboS"^) 
labeling the coherent state of the gravitational part of WDW theory, which is defined by 

dv e 5^e-*<'(''-^°)]w). (2.10) 

-CSO 

It also has the analogous properties of the coherent state of matter part, especially the resolution of the identity: 



The whole coherent state of WDW theory reads |^'2„)|^'r)o) = l^zo) ® \^vo)- On the other hand, due to the polymer- 
like structure, the coherent state of LQC is different from that of WDW. Here one can define Co = ;^^('^o + ibod^) to 
label the generalized coherent state [21, 34]: 



(*cJ :=X]e"^''^e-*''(''-''°)(t;|, (2.12) 
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where d is the characteristic width of the wave packet and 1 <C <C ?^ because of the semiclassical feature. For 
practical use, one defines the projection of this state on some lattice of variable v, saying the shadow state [34]: 



^yh^d._ J2 e-^^^^e^o^^^-^o'lfc), fceZ, (2.13) 
fe= 



where we chose the regular lattice {w = fc, fc G Z}. This shadow state also has the analogous properties of a coherent 
state. Note that our final result will not depend on the particular choice of regular lattice, since the "Hamiltonian 
operator" in Eq. (2.7) is a difference operator with step of size "4". The resolution of identity now reads 

where the identity I is in the subspace in which the states ^ have support only on the regular lattice. It should be 
noticed that the states with support on semilattice only in a single semi-axis of the real line was studied in Ref. 
[14], which are superselected by an alternative Hamiltonian constraint operator suitable to deal with some delicate 
issues in LQC. With this Hamiltonian operator of special symmetrized ordering, the singularity decouples in the 
kinematical Hilbert space and hence can be removed. It would be interesting to study also the effective equations 
of this Hamiltonian operator by coherent states supported on semilattice. However, we remark that, for practical 
calculations in our coherent state functional integral, the dynamical difference equation in [14] still needs to be suitably 
"simplified". 



III. COHERENT STATES FUNCTIONAL INTEGRALS 

In the path integral of the conventional non-relativistic quantum mechanics, one needs to compute the matrix 
element of the evolution operator e~^^*^ within the time interval At. However, the situation of cosmology of GR is 
very different since both WDW cosmology and LQC arc totally constrained systems, and the operator C is not a true 
Hamiltonian. Instead, we start from the physical inner product, i.e., the transition amplitude, of coherent states with 
normalization: 

To calculate the transition amplitude, we split a fictitious time interval At = 1 into N pieces e = and thus get 
gjaC _ giE"=i _ ]^^^^ e'^"^. Inserting N times of coherent states resolution of identity of \^zo) and Eq. (2.11) 
(or Eq. (2.14)), Eq. (3.1) can be casted into 

/oo 
da A--^*(vl/,^,f,JA|-(f,^,M',J, (3.2) 
-oo 

where 

N 



- r ds d^ r ^^^^ rr ^3 

- y_ d<t>.-^ --d^^ y_ • • • ^ 11 ||*.J|||*.„_J| ' (3.3a) 



with zn = Zf,zo = Zi,riN = ??/, and rjo = rji. Notice that the characteristic widths a and 5 at different steps are not 
necessarily the same. So wc have to denote (t„ and 5n in the semiclassical states and |5'r;„) respectively at the 

"n-step" . Now the main task is to calculate the matrix element of the exponential operators on coherent states. The 
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exponential operator e""'^ can be expanded as 1 + ieaC + 0{e^). For the purpose of a concise writing, we introduce 
some intermediate-step notations, 



P4>n = 



< + <-! " < + <-! 

Through a detailed calculation shown in the appendix, we get for the matter part 



N 



(vl/.Je-^|vl/,,_J 
i\ ll*.JIII*.„-JI 



where the inner product of two adjacent states is 



and hence the product of series (^'z„|'S'z„_i) can be expressed as 



K-i) 



(3.4) 



(3.5) 



N 



• exp 



Zn-l 



2(^^+1+^') 



' N 

n 



/ 2cr„i7„_i 



, (pn - (pn-l 

+ ^p^n — : — 



"4?l2 



(3.6) 



Here we introduced a virtual width aN+i by hand, satisfying ctat+i — cn = cn — fjv-i, in order to get the tidy sum in 
the exponential position. In the limit of A'' — >■ oo, cr jv+i will approach ajv = erf and hence does not effect the quantum 
dynamics. 

For the gravitational part, a careful calculation in WDW quantum cosmology shown in the appendix yields 



^ /iTr I pieaQwDW I iTr \ 

n \^rin\^ I ^Vn-ll 
l|VI> llll^ II 



AT 



n=l 



'Tin MM ^'rjn-l I 



iea ■ UttG 



N 



where w„ = '^""^2+^2""""' , &n = '^"'''g2^+jS '^n = ^4^17' ^'^'^ nn=i ||^^J|||^^^"„"_'j ^^'^^^ similar to 

Eq. (3.6). Now we take the limit N ^ 00 and substitute dr for X^^^^ e to get the functional integral formalism of 
the amplitude: 



A([%][*i]) =eKl^/l'-l^^l'+l''/l'-l''*l') J da J[V(P{T)][Vp4,{T)][Vv{T)][Vb{T)]e 



where 



dr i 



dr V4ct2 J ' dr \ Ah"^ ) h 
d [5%'' 



,+ 



-hv- aUnG 



(3.8) 



(3.9) 



(3.10) 



Here the "dots" over (p and v stand for the time derivative with respect to the fictitious time r. The functional 
measures are defined on continuous paths by taking the limit of TV — >■ 00: 



N-l 

n 

n=l 



27rS 



dVndbn 



(3.11a) 
(3.11b) 
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Ignoring the total derivatives with respect to r in Eqs. (3.9) and (3.10), we can read out the total effective Hamiltonian 
constraint in WDW quantum cosmology as: 



(3.12) 



Note that ^ and are the square of fluctuations of v and h respectively. They can be seen as quantum corrections 
to the leading term: v^h'^ — ivb. 

On the other hand, careful calculations in LQC shown also in the appendix give 



N 



AT 



exp 



- ieaSTTG^ ((tJ„2 + ^) (sin2(26„)( 



8 . 4 N 

dl + dlJ^ dl + dl_J 

2dl 



-iS\-D.{Ahn)Vn 



■(1- 



dl + dl_,^ dl + dl_. 



where Vn 



(3.13) 



Eq. (3.6). It should be noted that, in the above calculation, the inner product of two shadow states is 



'^Irfa ^ , and n^=i 111 ""'iiii^" also takes the form similar to 

n+^n-l II — Cnllll — <n-l" 



k' k 



■ — 5^72 \-tbn-l{k-Vn-l) 



\k) 



'"P 2{dl + dl_,) 
Using the so-called Poisson re-sum formula: 



+ ibn{Vn -v„-i) 



^exp 



S^-^Jll i{bn - bn-l){k - Vn) 

di 



k=—oo 

we have the summation about v as 



^ g{k + x)= e'"""" / dy ff(2/)e-^2-'=^ e Z, 

I J—oo 



k=—oo 



^exp (^-{k - Vnf/di - i{bn - ?>«-l)(fc - Vn, 

=?-/-.-'«■-- - 

k 



(3.14) 



and hence the inner product as 



{Vn - Vn-lf _ {bn - bn-l)'^dl 

'2{dl + dl_,) 4 



+ ibn{Vn -Vn-l) 



(3.15) 



where we kept only the fc = term, because in the continuous limit A'' — > oo, one has &„ — 6„_i — > and the terms 
corresponding to non-zero integer k are of tlie same or higher orders of O ie^^ 



and hence negligible under the 

semiclassical condition d^ 1. The matrix element (*.^^ |e*'^"®^<2° |^^^_^) can be obtained by similar method. Thus 
we obtain: 



^([*/][*i]) = e^l'^l'-l^'I'+l^^l'-l^-l') j da j[V(t>{T)][Vp4,{T)][Vv{T)][Vb{T)\t 



(3.16) 
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The effective Hamiltonian constraint in LQC can be read out as: 

(„2^^)(sin2(2&)(l-l) + |)+z^sin(46)(l-l) , (3.17) 



p1 1 



wherein the terms ^ and ^ are the square of fluctuations of v and sin(26) respectively. They are also quantum 

corrections to the leading term. 

At the first sight, both J^/fcff and J^^ g look problematic due to the imaginary part. One might even suspect the 
validity of the coherent state path integral in the models. However, a careful observation reveals that the real and 
imaginary parts of the leading terms can be synthesized into a Moyal-* product [31] in both models, i.e., 

vH"^ - ivb = ve'^i^'^''^'^") (bvb) =: v * (bvb), (3.18a) 
sin2(26) + iv sin (46) = ue* (^"^"-^^O ( gin (2b)v sin (26)) = ^; * ( sin {2b)v sin (26)) . (3.18b) 
Therefore the effective Hamiltonian constraint in WDW theory takes the form: 

^e. = -§ - ^ + 12.G (. * (6.6) + ^ + ^ + ^) , (3.19) 

while that in LQC becomes 

i^ff = - § - ^ + * ( sin (26). sin (26)) {l - ±) + '-^^^ " 1) + ¥ + ^ ' ^'-''^ 

To understand how the Moyal ^-product emerges in the gravitational part of the Hamiltonian, recall that both 

©WDW oc v{bvb) and 6lqc oc v(sin (26)'()sin (26)) are non-symmetric operators which can be regarded as a product 
of two self-adjoint operators. Thus, the coherent state functional integrals suggest the Moyal ^-product to express 
the effective Hamiltonian for the quantum system with a non-symmetric Hamiltonian operator. Now we explore 
the motivation of a non-symmetric gravitational Hamiltonian operator through the LQC prescription. It should be 
noted that the initial Hamiltonian constraint operator in LQC is actually self-adjoint in the kinematical Hilbert space 
[10, 24]. To resolve the constraint equation and find physical states, one feasible method is to rebuild the constraint 
equation as a Klein- Gordon-like equation and treat the scalar (f> as an internal time. As a result, the constrained 
quantum system was recast into an unconstrained system of non-relativistic particle whose dynamics is govern by a 
Klein- Gordon-like equation with an em,ergent time variable [10, 24]. The price to get this Klcin-Gordon-likc equation 
is that the new gravitational Hamiltonian operator Slqc becomes a multiplication of two self-adjoint operators, and 
hence it is no longer symmetric. But this does not indicate that one could not employ 6 in the intermediate step 
to find physical states. On the other hand, this non-symmetric Olqc just provides a suitable arena to examine 
the Moyal ^-product from the path integral perspective. The appearance of Moyal *-product in our path integral 
formalism indicates a possible duality between the path integral formulation on a non-commutative Moyal plane (see 
e.g. Ref.[35]) and the canonical quantization on an usual phase space with a "non- symmetric" Hamiltonian operator. 

We can also take another practical way to symmetrize Olqc at the beginning. For example, one can define a 
symmetric version of Glqc by 

©Lqc := ^(©LQC + e^Qc) « [«(shr(26)«shr(26)) + {siM^)vsiM^))v], (3.21) 

and then carry out the same procedure of above coherent state functional integral. In the calculation of matrix 
element (^^„|0LQcl^Cn-i)' ^® could think that the operators v and sin (26)wsin (26) in 6lqc act on bra \ and 

ket respectively, while sin (26)t)sin (26) and v in Q\^qq act on bra (^[.^^ | and ket |5[.^^_j) respectively. Then it is 

not difficult to see that the imaginary parts generated by Olqc and O^qq cancel each other. Hence for the symmetric 

Hamiltonian operator corresponding to ©lqCj we can get a real effective Hamiltonian constraint without imaginary 
part for LQC. Similarly, we can also get a real effective Hamiltonian constraint in WDW quantum cosmology by 
employing a symmetric version of Owdw- 
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IV. ON THE EFFECTIVE DYNAMICS 



Using the effective Hamiltonian constraints J^fj and which contain Moyal * -product, one may derive the 

corresponding dynamical equations. For LQC, we can define the evolution equations by: 



for any dynamical quantity f{v,b). Especially, the evolution of basic variables can be obtained as: 



(4.1) 



12ttG 



V = 



V* (i;sin(26)cos(26)(l-4e2)) + 



+ 



sin (26) cos(26)(l -4e2) 
2i2 

sin^ (25) (1 - 4e2) 



8£4 



2v{l - it) sin {2b} j * sin (26) + Avef 

-^(^^^(1 - 4£^) + {v' + ^)4sin^ (26)s^ - 2v\')d.e' 



(4.2a) 



(4.2b) 



where e = 1/d denote the quantum fluctuation of sin6, dbS^ = d{e^)/db and = d/dv. Similarly, we can use the 
effective Hamiltonian (3.19) in WDW quantum cosmology to get the evolution of basic variables as: 



b = 



_ 1277(7 

12nG 



6 62 \ 



2v*{vb) + -^+[-- — \d,s 



1 f9 v^e^ 



(2vb)*b^VE^--[-- 



(4.3a) 
(4.3b) 



where e = 1/(5 denotes the quantum fluctuation of 6. However, there seem no way to understand Eqs. (4.2a)- 
(4.3b) directly as effective classical equations because of the ^-product therein. To get physically predictable effective 
equations, we have to appeal to other possibilities. 

Since the Moyal * -product originates from the non-commutativity of operators, one can symmetrize the operator 
©LQC as Eq. (3.21) and repeat the procedure of coherent state functional integrals in last section. Then it is not 
difficult to get the effective Hamiltonian constraint for LQC as: 



- -i - + ^'K-^ 4) - + 1) • 



(4.4) 



which takes the same form as Eq.(3.20) but without the ^-product. Note that this effective Hamiltonian constraint is 
different from that obtained in Ref.[21] where a different Hamiltonian constraint operator was employed. Using the 
conventional Poisson bracket, we can get the evolution of v as 



V 



12ttG 



1 



« + ^ ) sin (26) cos (26)(1 - As') + [ - - [v' + ^] sin^ (26) 



1 



2£2 



sin^ (26)(l-4£2) 
8? 



.(4.5) 



Then a modified Priedmann equation can be derived as 

2 



1 + 



2w2e2 



sin (26) cos(26)(l-4£2) + 



- 1 + 



2w2£2 



sin^ (26) - 



sin^ (26)(l-4£2) 



H 2 



(4.6) 



where pc = g^^a^^^a is a constant. To annihilate sin (26) and cos (26) in Eq. (4.6), we use the constraint equation 
(4.4) to get 



sin^ (26) 



l-4£2 \p^E 



- 2e^ 



(4.7) 
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where p = ^5 is the density of matter, F = 1 + j^-r , and E = 1+ ^Ji^i ■ However, Eq. (4.6) looks problematic since 

it depends on the volume v of the chosen fiducial cell. This originates from the fact that we have to use the coherent 
states peaked on the phase points {v, b) in the path integral. In the final picture we have to remove the infrared 
regulator by letting the cell occupy full spatial manifold. In this limit, the irrelevant correction terms proportional to 
l/(t;e)^ could be neglected, while the relevant terms proportional to would be kept, since e was understood as the 
fluctuation of sin b which does not depend on the fiducial cell. We finally get 



2 _ SirGpc 



^LQC 




(4.8) 



where the positive and negative signs correspond to the expanding and contracting universe respectively. Since ^ is 
the square of fluctuation of p^, one has ^ 2 s <C 1. If one ignored the higher order quantum corrections, the modified 

Friedmann equation (4.8) would be simplified to the well-known form: i?LQC — ^'^^P — -^-^ However, Eq. (4.8) 
implies significant departure from classical GR, which is manifested in the bounce or re-collapse points determined 
by i^LQC = 0. For a contracting universe, the quantum bounce happens when 



1 - 4e2 Vpc ~ J 2 2 y (1 - 4£2)2 + {di,e^)2 ' 

Because e <C 1 and dbS^ ^ 1, the right hand side of Eq. (4.9) could be very close to but no bigger than 1. Hence the 
so-called quantum bounce of LQC will occur when p increases to pboun ~ Pc- On the other hand, for an expanding 
universe, the positive sign should be chosen in Eq. (4.8). As a result, the Hubble parameter would always keep 
non-zero unless approaches asymptotically. Assuming this is the case, in the asymptotic regime we would get 



^^^QC = ^ f - 2,^) (1 - As^ - (^F - 2,^ ) ) . (4.10) 



Therefore, under above assumption, a rc-coUapsc would occur if p decrease to pcoU ~ "^s^Pc which coincides with the 
result in the canonical theory [11, 21]. As pointed out in Ref. [21], in this case the inferred re-collapse is almost in all 
probability as viewed from the parameter space characterizing the quantum fluctuation s. 

For WDW quantum cosmology, one can also symmetrize the operator 6wdw and repeat the procedure of coherent 
state functional integrals. Then it is not difiicult to get the efltective Hamiltonian constraint: 

i-i+i^-°("'4)(*'+i)' 

which takes the same form as Eq.(3.19) but without the *-product. Using this Hamiltonian constraint, we can get 
the modified Friedmann equation for WDW cosmology as: 



2 SirGpc 



'WDW — 



(4.12) 



where the positive and negative signs correspond to the expanding and contracting universe respectively. It is obvious 
from Eq. (4.12) that there would be no bounce for a contracting universe. For an expanding universe, the Hubble 
parameter i?-^DW might vanish only if dbS^ approaches asymptotically. In this case, Eq. (4.12) implies that a 
re-collapse would also happen if the density of matter could decrease to p i=a PcoU = '^s'^Pc- Hence, once higher- 
order quantum corrections are included, the inferred re-collapse is a common effect in both WDW cosmology and 
LQC under the condition that quantum fluctuations approach constant asymptotically. Intuitively, as the universe 
expands unboundedly, the matter density would become so tiny that its effect could be comparable to that of quantum 
fluctuations of the spacetime geometry. Then the Hamiltonian constraint may force the universe to contract back. 



V. CONCLUDING REMARKS 



The minisuperspace models of quantum cosmology provide the good avena for testing the ideas and constructions 
of quantum gravity theories. A few physically significant results have been obtained in both WDW cosmology 
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and LQC. In LQC, the big bang singularity is resolved by the quantum bounce, and the effective Hamiltonian 
constraint with higher-order quantum corrections could even be obtained by the semiclassical analysis, which implied 
a possible effect of quantum gravity on large scale cosmology. It is desirable to study such kind of predictions from 
different perspectives and in different frameworks. Since the higher-order corrections of the Hamiltonian come from 
the quantum fluctuations, a natural attempt to achieve them is to employ coherent state path-integral. On the 
other hand, the so-called Moyal * -product in quantum mechanics is generalized to many situations including quantum 
spacetime itself. It is also possible and desirable to derive the * -product by coherent state functional integral approach 
within quantum cosmological models. These issues have been addressed in previous sections. 

We summarize our main results with a few remarks. First, by the well-established canonical theories, the coherent 
state functional integrals for both WDW cosmology and LQC have been formulated by group averaging. As far as 
we know, this is the first attempt to apply coherent state functional integral to the models of quantum cosmology. 
Second, the main calculation results of our coherent state functional integrals are the effective Hamiltonian constraints 
(3.19) and (3.20) for WDW cosmology and LQC respectively. These show that the Moyal (star) -product can emerge 
naturally in the path integral approach via the effective Hamiltonian with higher-order quantum corrections. That 
is the main reason why we start with the non-symmetric gravitational Hamiltonian constraint operators ©wDW and 
©LQC for the path integrals. Whether the resulted Hamiltonian constraints with *-product could make some physical 
prediction would be an interesting open issue. Tentatively, the appearance of Moyal ^-product in our coherent state 
path integral indicates a possible duality between the path integral formulation on a non-commutative Moyal plane 
and the canonical quantization on an usual phase space with a "non-symmetric" Hamiltonian operator. Third, for 
symmetric Hamiltonian constraint opcirators, the effective theories and modified Friedmann equations have been 
obtained by the coherent state path integrals in both WDW cosmology and LQC. For LQC, the effective equation 
(4.8) can reduce to the first-order modified Friedmann equation when higher order quantum corrections are neglected. 
Hence the quantum bounce resolution of big bang singularity can also be obtained by the path integral representation. 
On the other hand, if higher order corrections are included, under the condition that quantum fluctuations approach 
constant asymptotically, there is great possibility for the re-collapse of an expanding universe due to the quantum 
gravity effect, which coincides with the result obtained in canonical LQC. Moreover, the effective equations imply 
that the inferred effect of re-collapse is common in both WDW cosmology and LQC under above condition. Finally, 
it should be noted that, as we used the coherent states of Gaussian type, the effective equations and hence the 
inferred effect of re-collapse are only valid with the assumption that these coherent states can faithfully represent the 
semiclassical behaviors of WDW cosmology and LQC. Whether there is a similar result for other semiclassical states 
is still an interesting open issue. 
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Appendix A: Calculation of the functional integrals 



We will use the following Gaussian integrals: 

r+oo 



f 

J — c 



dx e 



+00 



a 2a2' 



dx cos(bx)e " ^ = e 
a 



/ 
/. 

/+00 
dx x^ sin {bx)e~ 
-00 

/+00 
dx x^ cos {bx)e~ 
-00 



dx X cos (6x)e " ^ = — e 

aa; a;sm (oa;)e = 4^, 



a V4o^ 



63 \ 


b2 




e i^, 


362 


64 \ 


4a6 " 


^ 16a8 i 



e ^o'i 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 
(Alf) 
(Alg) 
(Alh) 



with constant parameters a > 0, 6 € M. For the simplicity of notation, we denote the momentum conjugate to ^ by p 
rather than in this appendix. 

The main task is to calculate the matrix elements of exponentiated operators: (\l/2^|e*'^"R^|\l/2„_i), 
(*^„|e-'^"®w°w|^^^_j and (*.^„|e-'^"®i'«'^|*^^_ J . For matter part, we get 



AT, l,T, \ / y-Pn - (l^n-lf ^ i CrlPn + <7l-iPn-l r , , x 

(*^„|*^„_i) = exp ( , _^ ^ + _2 , _2 (</'n - 



>exp 



Introducing some intermediate-step notations: 



- T{Pn~ Pn- 



-7 ^ gn-i</'n 



we can get a concise writing: 

{4>n - 4>n-lY , i 



Pn — _9 , _2 ' ''^n — _9 , _2 



(*^„|*^„_,) =exp 



)exp 



— " ] exp \^--^{pn-Pn-l){<f>- ' 



2{al + ol_^) ' h 

and the integration wherein is easy to be done by changing the integrated variable (j) to <p = (j) — (j>^: 

(Pn -Pn-l)^0-2 \ 



'COS I -^{Pn -Pn-l) 



'Try cr„ exp 



4?l2 



Finally we get the inner product: 



t>n-l) (Pn- Pn-l) 0-2 i 



2(^^ + <-i) 



4?i2 



+ J^Pn{<i>n - <i>n-l) 



(A2) 
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To get Eq. (3.6), we need to deal with the product 



2<Jn(T„_i 

2 _L^2 I exp 



TV 



■ exp 



E 

U=l 



.n=l 



where the summation in the last part of exponential in Eq. (A3) can be re-organized as 



(A3) 



N 

E 



'Pi 



0Af - 2(pM<pN- 



2{al + al) 2{al + al)J 2{al + al) 



N 



^ 2(0"n+l + (^n)M<f>n - (f>n-l) - (^n+l + Cr„_i)((7„+i - (Tn-l)(f>l ^^^^ 



n=l 



and similarly 



JV 

E 

n=l 



{pI - 2pnPn-lWl pI-i(^1 



Pn'^n+i^n 



2 2 2 



4^2 



2hHa%^,+a%) 2^2(^2+^2) 



AT 

E 

n=l 



4('7n+lO'nO'n-l + Cr^P-2_^ )p„ (p„ - p„^i) + 2cr^(cr„+i + i7„_i)(o-„+i - CTn-Qp^ 

4S2(a2 + i+^2)(^2 +^2_^) 



(A5) 



Collecting the above results of Eq. (A3), Eq. (A4) and Eq. (A5), we can finally get Eq. (3.6). The matrix element 

-2 

(^'z„|zea|7|\E'z„_i) is proportional to 



(*^„|P^I*^„-i) =exp 



■I 



dcj) 



■ 2 fi" 

Pn-l + ~T 



>n-l "n-1 



2(a2+a2_i) ■ h 

(A A. ^2 , 2ihpn-l ^^ . 
"4 [(p-9n-l) H 2 {(p-9n-l) 



exp 



-=^-;^(Pn-Pn-l)(<^-<^J )• 



To do the above integral, we have to rewrite the integrand as a function of </> — Except for the exponential function, 
the left terms are polynomials of </>. The integral of zeroth power term p'^_i + ^ gives 



V^V exp ( - 



(p„ -Pn-l)2g2 
4/12 



• Pn-l + 



and hence its contribution to matrix element (^z„ |p^|*z„_i) is: ('J'z„|^'z„_i) (pn-i + To deal with the linear 

and square terms of (j), we first rewrite 
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and get the integrals: 



#(0 - 4>n? exp - ^'^ 4"^ - -(pn - Pn-lM - 4>n) 



■v/tt Y cr^ exp 



{Pn - p„-l)V2 W 0-2 (p„ - p„_i)2(cr2)2 



4^12 



J #(^- '?^„)exp 1^- ^'^ - -p„_i)((?i- 0„) 



-v/tt V 0-2 exp 



CT, 

{Pn-Pn-1 



f<^n\ -APn-Pn-l)<yl 



2h 



Combining all the above results, the matrix element znW^°'~^ \^ zn-i) 

;i2 ;i2^ (p„-p„_i)2(^)2 ;i2(^^-)2 



zea , , , / o ft 



ih(j)„al{Pn - Pn-l) Pn-K^liPn - Pn-l) 2ihpn-l(, 



'n-1 



n-1 



(*^J*^„_i)exp 



( Pn-1 + 



^2 



1 pmatt 



(A6) 



up to C'(e^)- Here denotes a polynomial of pn — Pn-i and 0„ — 4>n-i without zcroth order terms. To get a 

functional integral formalism, one has to take the number of steps TV — > oo and so (pn — 4'n-i — > 0, Pn—Pn-i — > 0. 
Furthermore, there is already an infinitesimal e multiplied to P™*^**. This means that P^n-i does not become a time 
derivative in the limit of N — > oo. As a result, P^n-i does not contribute to effective action 5™^". 
For gravitational part of WDW theory, we get 

(*^Je-^^«®wDw|^^^_^) = (vl/^jM/^„_,) +iea(127rG) j dv%^{v)vd,[vd,^ r,„_,{v)\ +0{e^). (A7) 
Similar to the matter part, it is easy to get the inner product of two coherent states as 

/,T, I.T, \ /- fjT ( {Vn-Vn-lf (6„ - 6„_i)2(52 ^ 



(A8) 



where some notations are defined as before: 



The second term of Eq. (A7) is proportional to 



-^^;^S^+i(6„-6n-l)K-^7n) 

e 



"n-1 



- ihn-\v 



(A9) 
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where w„ = ^" ^2Ttf"'"" ^ • To do the integral, we first have to rewrite the integrand as a function oiv — Vn as follows: 



1 , 1,2 \ 2 , v'^{v-Vn-lY , '^ihn-iv'^ {V - Vn-l) v{v - Vn-l) 



+ K-1 + ibn-lVn 



^1 ^ ^l-l 

{V - VnY 



51-^ 



+ 



(u — t; 


n- 




1 




L 






1 




"n-1. 


( 



- ihn-iv 



2Vr, 



n .2 



{V - Vn) 



"n-1 



5U 



+ ibn- 



(v - Vn) 



+ 



+ 



where = Vn-i = 5i+s'i_^ — 



j:4 A2 S"^ I 

"n-1 "n-1 "n-l/ 

Doing the integral term by term, we can obtain the result of Eq. (A9): 



(AlO) 



j dv^*Jv)vdy[vdy^nn-iiv)] = -(*r,„|*r,„_i) (^-^ + bl_-^^vl + ibn-lVn 



+((W-+Ci)-^^-%^+W-)f (All) 



and finally the matrix element: 

2« 



i6„_,I;„(|^-l)+P™)j^ (A12) 



where -P^,^ is a polynomial of ii^ — i^n-i, 6„ — and J„ — (5„_i without the zeroth order term. As before, this 
quantity does not contribute to the effective action under the continuous limit e = 1/A^ — \ 0. 

For gravitational part of LQC, we have to calculate the inner product (4!.^^ and matrix element 

iea0LQc|*.,^„_i)- '^'^^ inner product of two shadow states has been obtained in Eq. (3.15). The order of 

0{e) of the matrix element (*f^|e"""®^«<^|*^^^) is 

(*^J-ieaeLQc|*Cn-i) = iea^E [^(^ + 2)*CnW*C„-i(^ + 4)-2fc'*^„mc„_,(fc) 

k 

+k{k - 2)**^ (fc)*c„_i {k - 4) 
^ iea^ - Dl^_, + D-^_,) . (A13) 

Now we need to deal with the three terms n-n ^n.n-ii ^n.n-i separately. First, wc get 



£'^n-i = E(^' + 2fc)e 



(fc-^„)^ (fc + 4-^„_l) 
2dH 2dH 



= exp 



4^{Vn - Vn-l) 



+ i^bn - 



1 ^ — ibn{k — Vn)+ibn-l{k+'i—Vn-l) 
{Vn - Vn-l)"^ 



+ 2fc) exp (-(A - ZJ+) V4 - t{bn - bn-i){k - U+)) , 



+ ibn{v n Vn— 1 



(A14) 



where v^ = ^^^2^^^!^'" ^ To do the summation in the above equation, we first have to rewrite + 2fc as a 
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function of k — vt: 



e + 2k = (fc-tJ+)2 + (2tJ^ + 2)(fc-^J+)+^jf + 2tJ+ 



= {k-v+y + 2{u 



+ 



where Vn = ^"^^rfa""" ^ • The sum about the square term is 



^(A: - v^f exp (-(fc - ^7+)Vci2 - ^{hn - 6„_i)(fc - ^7+)) 

k 

= V^Y^cxp (-(6„ - bn-ifWjA) idi/2 - ibn - 6„_i)2(d2;)2/4^ . 
The sum about the linear term is proportional to 

Y,{k - vt) exp (-{k - v+f/d^ - i{hn - bn-l){k - vl)) 

k 

= -iVT^\/(^exp (^-{bn - bn-i fdl/4j [bn - bn-i)dl/2, 
and the contributions of the zeroth order terms are proportional to 

J2 exp {-{k - v^fK - i{bn - bn-i){k - vi)) = yiv^exp (-(6„ - bn-lfWj^ . 
k 

Collecting the above results, we get 

2dl 



(Vn) - 2Vn 



dl + dl_, 2 



dl^p^ 

n,n-l I ' 



(A15) 



where P^n-i denotes a polynomial of Vn — Vn-i, bn — bn-i and d„ — d„_i without the zeroth order term. Here we 

have expanded the factor exp (^"^^jq^^jr— in Eq. (A14) as 1 — ^^Ja"^^;^""^^ H • Except for the leading term 1, all 

the other terms can be conflated with P^ n-i- Analogous to the matter part, under the continuous limit A'' — > oo, 
this Pn.n-i does not contribute to the effective action of gravity. With the experience of computing Dt^_-^, it is easy 
to calculate -D° ^-i ^^d D~^_-^ as follows: 



dl 



^dl ,d^,p- \ 

dl + dl_,^ 2 ^ "'"-M- 



(A16) 
(A17) 



Taking the expansion e 
get the combination 



V = 1 _ ^ (^) neglecting the higher order terms than (^), we can 



n+ — nO -\- 



= -4(*cJ*a-i) [((^")' + f ) (si^' (2S»)(1 - dj+fc) + dj+iir) + »sin(46„K(l 



and hence the matrix element {^_i^J\e~^'^°'®^'^'^\^^_^) is 

/ - S 4 

.»„ . 3.g(((.„)» + ^){*' (25„)(1 - -^j^) + -^^^ 

« 2dl 



P, 



n,n—l 



(*cJ*C„_i)exp 



+i sin (46„)t;„(l 



I pgrav \ 

dl + dl_,'dl + dl_,^ "'"-V 



(A18) 
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